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We have derived the closed system of averaged MHD-equations for general oscillat- 
ing flows, which are purely oscillating in the main approximation. We have used the 
mathematical approach which combines the two-timing method and the notion of the 
distinguished limit. Properties of the commutators are used to simplify calculations. The 
direct connection with a vortex dynamo (or the Langmuir circulations) has been demon- 
strated and a conjecture on the MHD-dynamo has been formulated. 



1. Introduction 

This paper derives the averaged MHD-equations for oscillating flows. The resulting 
equations are similar to the original MHD-equations, but surprisingly (instead of com- 
monly expected Reynolds stresses) the drift velocity (or just the drift) plays a part of an 
additional advection velocity. 

It is known that the drift can appear from either Lagrangian or Eulerian consid- 
erations. The Lagrangian drift appears as the average motion of Lagrangian particles 
and its theory is often based on the averaging of ODEs, see Stokes (1847)| Lamb (1932) 



Longuet-Higgins (1953)||Batchelor (1967)[|Andrews fc Mclntyre (1978)[|Craik (1982)||Yudovich (2006) 
In this paper we focus on the Eulerian drift, which appears as the results of the Eulerian 
averaging of related PDEs without addressing the motion of particles, see Craik fc Leibovich(1976)| 
Craik (1985) [ |Riley (2001)] | Vladimirov (2010)||IlTn~&; Morgulis (20lT)| The detailed ma- 



terials about the Eulerian drift can be found in Vladimirov (2010) 



To derive the averaged equations we employ the two-timing method, see e.g. Nayfeh (1973) 



Kevorkian & Cole (1996) We expose it as an elementary, systematic, and justifiable pro- 
cedure that follows the form developed by Yudovich (2006)| Vladimirov (2005) Vladimirov (2008) 



Vladimirov (2010) This mathematical procedure is complemented by a novel material 



on the distinguished limit, which allows to find the proper slow time-scale. 



2. Functions and operations 

We introduce functions of variables x = (xi, X2, £3), s, and r, which in the text below 
serve as dimensionless cartesian coordinates, slow time, and fast time. 
Definition 1. The class H of hat-functions is defined as 

/eH: /( !B) 8,r) = /(aj J s J T + 27r) (2.1) 

where the r-dependence is always 27r-periodic; the dependencies on x and t are not 
specified. 
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Definition 2. For an arbitrary / € H the averaging operation is 



1 



T0 + 27T 



(f) = ^J f(*,*,T)dT, Vr (2.2) 

where during the r-integration s = const and (/) does not depend on tq. 
Definition 3. The class T of tilde-functions is such that 

/GT: f(x,s,T)=f(x,s,T + 2ir), with (/) = 0, (2.3) 

The tilde-functions are also called purely oscillating functions (T-function represents a 
special case of H- function with zero average) . 

Definition 4- The class B of bar-functions is defined as 

/GB: / T =0, J(x,s) = (f(x,s)) (2.4) 

(any H- function can be uniquely separated into its B- and T- parts with the use of (|2.2p ) 
Definition 5. T -integration (or tilde-integration): for a given / we introduce a new 
function f T called the T-integral of /: 

r = J f( x > s > a ) da -^l f(x,s,a)do)dn (2.5) 

which represents the unique solution of a PDE df T jdr = f (with an unknown function 
f T and a known function /) supplemented by the condition (/ } = (f T ) = (|2.3p . 

The r-derivative of T-function always represents T-function. However the r-integration 
of T-function can produce an H-function. An example: / = f 1 smr where f 1 is an 
arbitrary B: one can see that (/} = 0, however (JJ" f(x,s,cr)dcr) = /j ^ 0, unless 
f 1 = 0. Formula (|2.5[) keeps the result of integration inside the T-class. T-intcgration is 
inverse to T-diffcrcntiation (f T ) T — (f T ) T — f\ the proof is omitted. 

Definition 6. A dimensionless function / = f(x,s,r) belongs to the class O(l) 

/ 6 O(l) (2.6) 

if / = 0(1) and all required partial x-, s-, and r-derivatives of / are also O(l). 

Here we emphasize that through all the text below all large or small parameters are 
represented by various degrees of a only; these parameters appear as explicit multipliers 
in all formulae containing tilde- and bar-functions; these functions always belong to O(l)- 
class. 

We also will use some properties of r-derivatives such as 

7r = /r, (/r) = (/r)=0 (2.7) 

The product of two T-functions / and g represents a H-function: fg = F, say. Separating 
T-part F from F we write 

F = F-(F)=fg- (fg) = {fg} (2.8) 

where the notation {•} for the tilde- part is introduced to avoid two levels of tildes. We 
will use that the unique solution of a PDE inside the tilde-class is 

df/dr^Q 7=0 (2.9) 

which follows from (|2.5p . Since the average operation (|2.2p is proportional to the inte- 
gration over r, then from the integration by parts we have 

<[a,b T ]) = -<[5r,6]) = -<[a T ,b]), ([5,M) = -([a T ,b]) = -([a T ,b]) (2.10) 
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where [a, b] stands for the commutator of two vector fields a and b which is antisymmetric 
and satisfies Jacobi's identity for vector fields a, b, and c: 

[a,b]=(b-V)a-(a-V)b, (2.11) 
[a,b]=-[b,a], [a,[b,c]] + [c,[a,b]] + [b,[c,a]]=0 (2.12) 

A useful property of the commutator is 

diva = 0, divb = => div [a, b] = (2.13) 

For any tilde- function a and bar- function b (|2.10l) . (|2.12[) give 

([a,[b,a T ]}) = [b,V a ] where V a = i([5,5 r ]} (2.14) 

3. Two-timing problem and distinguished limits 

The governing equation for MHD-dynamics of a homogeneous inviscid incompressible 
fluid with velocity field it*, magnetic fields h* , vorticity u>* = V* x it* and current 
j* = V* x h* is taken in the vorticity form 



[w*,it*]* - \j*,h*)* = 0, in P* (3.1) 
[h*,ur = 0, 



or 
dh 

dv 

V* • it* = 0, V* • /i* = 



where asterisks mark dimensional variables, £*-time, x* = (x\, x^-cartesian coordi- 
nates, V* = (d/dxl,d/dx2,d/dx^), and [•,•]* stands for the dimensional commutator 
(|2.1ip . In this paper we deal with the transformations of equations, so the form of flow 
domain T>* and particular boundary conditions can be specified at the later stages. 

We accept that the considered class of (unknown) oscillatory solutions u* , h* possesses 
characteristic scales of velocity U, magnetic field H , length L, and high frequency a* 

U, H, L, ct*>1/T; T = L/U (3.2) 

where T is a dependent time-scale. In the chosen system of units the dimensions of U and 
H coincide; we choose them being the same order U = H. The dimensionless variables 
and frequency are 

x = x*/L, t = t/T, u = u*/U, h = h*/U, a = a*T > 1 (3.3) 

We assume that the flow has its own intrinsic slow-time scale T s i ow (which can be different 
from T) and consider solutions of (|3.ip in the form of hat-functions (|2.1I) 



u* = Uu{x,s,t), h* =Uh(x 7 s 7 T); with r = at, s = fi t, fi = T/T slow (3.4) 
Then the use of the chain rule and transformation to dimensionless variables give 

^ + ^)Q + -[Q,u]--\j,h]=0 (3.5) 
OT a OS J a a 

d ft d \ t 1 rC _ n 
^T + -7T )h+-[h,u}=0 
or a os J a 

div h = 0, div u = 



In order to keep variable s 'slow' in comparison with r we have to accept that H/cr <C 1. 
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Then eqn. ([3.5)1 contains two independent small parameters: 



11 lfi , , 

la* a JsiowCT* a 

Here we must make an auxiliary (technically essential) assumption: after the use of the 
chain rule Q3.5)) variables s and r are (temporarily) considered to be mutually indepen- 
dent: 

t, s — independent variables (3-7) 

From the mathematical viewpoint the increasing of the number of independent variables 
in a PDE represents a very radical step, which leads to an entirely new PDE. This step 
should be justified a posteriori by the estimations of the error of the obtained solution 
(rewritten back to the original variable t) substituted to the original equation (|3.1|) . 

In a rigorous asymptotic procedure with a — > oo one has to consider asymptotic paths 
on the (e,£i)-plane such that 

(£,£!)-► (0,0) (3.8) 

Each such path can be prescribed by a particular function Q(a). One may expect 
that there are infinitely many different (although some of them can coincide) solutions 
to ([3.5)1 corresponding to different fi(cr). However for these equations (as well as for 
many others) a unique path can be found, which is called the distinguish limit (or the 
distinguished path). The notion of the distinguish limit is practical and heuristic, sec 
Nayfeh (1973)] |Kevorkian fc Cole (1996)] its definition can vary for different equations 



and in different books and papers. For our problem we write that the distinguished limit 
is given by such a function SI = Sl^cr) that allows to build a self- consistent asymptotic 
solution. Here the term self- consistent asymptotic solution means that the required suc- 
cessive approximations can be calculated. These calculations include the elimination of 
the reducible secular in s terms; the reducible secular terms are such terms which can be 
excluded by increasing the slow time-scale. For instance, a non-secular term proportional 
to sins gives secular terms in its Taylor's decomposition with respect to t = as). Below 
we show that the choice 

Cl(a) = l/a: r = at, s = t/a, a = const (3-9) 

allows to build the distinguished limit solution. The uniqueness of such a path for the 
considered class of solutions can be proven but we avoid such details in this paper. Hence 
the governing equations are 

u T +e[Q,u] -e\j,h] +e 2 Q s = 0, £ = l/a^0 (3.10) 
h T + e[h, u] + e 2 h s = 
div u = 0, div h = 

where the subscripts r and s denote the related partial derivatives. 



4. Derivation of the MHD-Drift averaged equation 

Let us look for solutions of p.lOp in the form of regular series 

oo 

(h,u) = ^ £ fc (^,wfc); h k ,u k emno(i), * = 0,1,2,... (4.1) 

k=0 
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In this paper we enforce the restriction 

u = 0, h Q = (4.2) 

which is natural physically if one considers, say, how the secondary vorticity develops on 
the background of a wave motion. The substitution of (|4.ip . (|4.2p into (|3.10p produces 
the equations of successive approximations. The equations of zero approximation are 

w r = w 0t = 0; h 0T = h aT = (4.3) 

Their unique solution (|2.9[) is u>o = and ho = 0. Taking into account (|4.2[) we can write 

£D = 0, h Q = (4.4) 

which means that in zero approximation the flow is potential, purely oscillating, and the 
magnetic field vanishes. This leads to the similar equations for the first approximation 
of (l3~T0l) . (l4T|) -(l4a 



Q>ir=0; h lT = (4.5) 

which have the unique solution 

U!=0, /m=0, wi=0, ^1=0 (4.6) 
where mean functions remain undetermined. The equations of second approximation that 



take into account (|4~2|) . (|4~4|) . (|4~6 |) are 

Z>2t + [Wi,«o] = 0, U>2r + [Wi,« ] = (4.7) 

which after T-intcgration (|2.5|) yield 

W 2 = K,Wl], h2 = [u T ,h 1 ], ^2=0, ^2=0 (4.8) 

The equations of third approximation that take into account (|4.2p . ([4.4p . (|4.6[) arc 

w 3t +u 1s + [w 2) «o] + - = (4.9) 

/^3r + hi s + [h 2 , Uq] + [h!,Ui] = 

The bar-part (|2.2p of this system is 

uJts + [wi.ui] - + ([w a ,«o]> = (4.10) 

ft,i s + [hi, ui] + ([h 2 ,2o]) = 
which can be transformed with the use of (|4.8p and (|2.14p to the final form 

u u + pi,ui + Vo}-tf 1 ,hi} = (4.11) 
h ls + {hi,ui + Vo\ =0 

Vo^\([u ,u T Q }) (4.12) 

If one uses these equations as a closed mathematical model, then all the subscripts and 
bars can be deleted: 

u> s +[u>,u + V]-\j,h]=0 (4.13) 
h s + [h, u + V) = 

V=±{fr,tT]) (4.14) 

One can see that: 
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1. The equation for the oscillating velocity uq in our consideration is absent, there 
are only two restrictions: Uo is incompressible and potential. Hence the drift velocity Vq 
(I4.12p represents a function that is 'external' to the equations. 

2. The derived system of equations (|4.11[) looks similar to the original one (|3.1|) . One 
may think that (|4.11|) describes 'just' an additional advection of vorticity and magnetic 
field by the drift. However, the fact that the averaged vorticity is transported with such 
an additional velocity is highly non-trivial; in particular, it contains the possibility of the 
vortex dynamo or Langmuir circulations, which we consider below. 



5. Stokes drift 

Our description of the drift differs from a classical one, therefore we first demonstrate 
the match of (|4.12|) with the classical Stokes drift. Let velocity field Uq and £ = Uq be 



Uq(x, s, r) = p(x, t) sinr + q(x, t) cost 
£ (x, t, t) = —p(x, t) cost + q(x, t) sinT 
with arbitrary B-functions p and q. Straightforward calculations yield 

[«o,£ ] = [P. 9] 

hence the commutator is surprisingly not oscillating. The drift velocity (|4.12[) is 

Vo = ^([«o,C ]) = 

The dimensional solution for a plane potential harmonic travelling wave is 

cos(fc*x* - t) 
sin(fc*a;* — r) 



(5.1) 

(5.2) 
(5.3) 

(5.4) 



Uq = UUq, , Ito = CXp(fc*Z*) 

where (x* , z*) are cartesian coordinates and k* = 1/L is a wavenumber. In Stokes (1847) 



Lamb (1932) Debnath (1994) one can see that U = k*g*a*/a* where a* and g* are 
dimensional spatial wave amplitude and gravity; however these physical details are ex- 
cessive for our analysis. The dimensionless velocity field (|5.4f) and £ are 



u 



cos(x — r) 
sin(x — t) 



- sin(x 
cos(a; - 



■T) 

T) 



(5.5) 



where both fields (|5.5[) are unbounded as z — > oo, but it is not essential for our purposes. 
Fields p(x,z), q(x,z) ()5.1j) are 

(5.6) 



(5.7) 



p = Ae 

The calculations (with the use of (|5.2[0 yield 



/ sin x \ 


, q = Ae z [ 


cosx \ 


V —cos a; / 




sin a: / 



Vr t = 



,2; 



The dimensional version of (|5.7|) is 



U 2 k* 



2k* z* 



(5.8) 



which coincides with the classical expression for the drift velocity given by Stokes (1847) 



Lamb (1932) Debnath (1994) To obtain (|5.8j) one should take into account the difference 



between the time t and s = t/a (|3.9 
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6. The averaged Euler's equations and vortex dynamo 

A special case of (|4. . (|4. 13[) without a magnetic field is 



V] = 0. 



divw = (6.1) 
Different s-independent versions of eqn. (|4.11j) were derived in the studies of Langmuir cir- 



culation by Craik & Leibovich(1976) and for the steady streaming problems by Riley (2001) 
Ilin fc Morgulis (2011)| the methods employed by these authors are different and more 
cumbersome than our method. In order to demonstrate the possibility of vortex dynamo 
we first notice that eqn. (|6.ip can be integrated (in space) as 

u s + (u ■ V)« + oj x V = -Vp, V • u = (6.2) 

where p is a function of integration and the second equation follows from the continuity 
equation in ()3.1j) . Let the zero approximation (|4.4|) represent the plane potential travel- 
ling gravity wave (|5.5p with the drift velocity (|5.7p . Let cartesian coordinates (x, y, z) be 
such that V = (U, 0,0), U = e 2z , u = (u,v,w) where all components are ^-independent 
(translationally-invariant) 
ponent form of (|6.2|) is 



and 



-variables coincide with ones in (15.51). Then the com 



UiUy 

vw„ 



wv. 



WW 



= 

- UUy - 
y — UU Z 



-Py 



which can be rewritten as (see Vladimirov (1985)] Vladimirov (1985a)) 



V S + VVy + WV Z = —Py — pQy (6-3) 

W s + VW y + WW Z = — P z — p<& z 

Vy + W z = 

P S + Up x + VPy = 

where p = u, $ = U = e 2z . and P is a modified pressure. One can see that (|6.3j) is 
mathematically equivalent to the system of equations for an incompressible stratified 
fluid, written in Boussinesq's approximation. The effective 'gravity field' g = — V$ = 
(0, 0, — 2e 2z ) is non-homogeneous that makes the analogy with a 'standard' stratified fluid 
non-complete. Nevertheless one can see that any increasing function u{z) = p(z) (taken 
from the shear flow (u,v,w) = («(z),0,0)) produces 'Taylor instability' of an inversely 
stratified equilibrium. It leads to the growth of longitudinal vortices and can be connected 
to Langmuir circulations, see Craik fc Leibovich(1976)| Leibovich (1983)| Craik (1985) 
Thorpe (20047 



7. Discussion 

1. The consideration of this paper is based on the assumption that the enforced fre- 
quency a* (|3.3[) of oscillations is higher than all intrinsic frequencies. This frequency 
appears in our theory via the prescribed potential velocity u (|4.3p . 

2. The prescribed oscillatory velocity Uq can be caused by different factors. For exam- 
ple, it can be produced by oscillations of boundaries or appear in full viscous theory after 
the matching of external flow with boundary-layer solution. The latter option is often 
considered, see |Riley (2001)| |Vladimirov (2008")] |Ilin fc Morgulis (2011) 

3. To justify the distinguished limit (|3.9p mathematically, one should prove that any 
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different path 0(cr) produces an asymptotic solution which contains terms secular in s or 
does not produce any asymptotic solution at all. The following statement can be proven: 
for Vq 7^ (|4.12[) and the function f2(<r) = l/tu a (with a constant a > —1) all solutions 
with a < 1 contain secular terms, while all equations with a > 1 produce a controversial 
(unsolvablc) equations of successive approximations. If Vq = then the statement is 
different but we do not describe it here. The proof is omitted. 

4. There is a challenging physical fact to explain: the existence of the distinguished 
limit (|3.9p means that there is a hidden slow time-scale T s i ow = a*T 2 in the system. 

5. The consideration of translationally-invariant MHD-motion in (|4. is possible in 



the spirit of the analogy between MHD flows and stratified flows, see Vladimirov, Moffatt and Ilin (1996) 



6. The mathematical justification of the equation (|4. 1 1[) by the estimation of the error 
in the original equation (|3.1|) is easily achievable. 

7. The higher approximations of the averaged equation ([4.11)1 can be derived. They 
are especially useful for the study of motions with Vq = 0. In particular, one can show 
that in this case Langmuir circulations cab be still generated by a similar mechanism. 

8. The viscosity and diffusivity can be routinely incorporated in ([4.11[1 as the RHS- 
terms z/V 2 aJi and nV 2 hi . Accordingly, viscous and diffusion terms will appear in the 
equations ([6.2)1 and (|6.3j) . At the same time after the incorporation of viscosity one 
more small parameter appears in the list (|3.6[) . and the distinguished limit should be 
reconsidered. 

9. The incorporation of the density stratification and gravity field into presented theory 
(or as a separate theory) is straightforward. 

10. The abolishing of the requirement of a vanishing mean flow in zeroth approximation 
(j4.2|) is also straightforward. However in this case the distinguished limit ()3.9|) is different 
and the resulted averaged equations are more complicated than ([4. lip . 

11. For the finite and time-dependent flow domain D{t) the definition of average ([2.2)1 
directly works only if x G V at any instant. If at some instant x ^ V then the theory 
should include a 'projection' of the boundary condition on the 'undisturbed' boundary. 
Such a consideration requires the smallness of the amplitude a* / L of spatial oscillations 
of fluid particles. However one can see that a* ~ u* /a* and hence a* /L ~ l/a = e ([3.6)1 . 
Therefore the considering of a time-dependent domain does not introduce any new small 
parameter, and the distinguished limit Q3.9P will stay the same. In particular, it means 
that our small parameter e is the same as the dimensional slope of free surface in the 



theory of Langmuir circulations by Craik & Leibovich(1976) 



12. The determining of the function u(z) = p(z) in ([6.3p (for real Langmuir circulations) 
requires an additional theory which includes viscosity or turbulent tangential stresses at 
the free surface. It is interesting, that from this viewpoint an 'unstable stratification' can 
be continuously generated and amplified by tangential stresses applied at free surface. 



We do not compete here with the theories by Craik & Leibovich(1976) Leibovich (1983) 
Thorpe (2004)[ which describe this phenomenon well 



13. One can suggest that since the equations ([4. lip for h = do describe a mecha- 
nism of vortex dynamo, and the mathematical structure of the full averaged equations 
(|4.1ip . ([4.13p with h ^ is similar, then these full equations could also describe a possible 
mechanism of MHD-dynamo, such as the generation of the magnetic field of the Earth. 
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